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1: Introduction 


3: Lidar variables dependency on refractive index and size 


5: Simulated lidar results based on dust model 


The physical and chemical properties of soil dust aerosol particles fundamentally affect their 
interaction with climate, including 

• Shortwave absorption and radiative forcing, 

• Nucleation of cloud droplets and ice crystals, 

• Heterogeneous formation of sulfates and nitrates on the surface of dust particles, 

• Atmospheric processing of iron into bioavailable forms that increase the productivity of marine 
phytoplankton. 

To distinguish between aerosols with different physical and chemical properties, lidar 
measurements are frequently used, such as 

• Extinction-to-backscatter (lidar) ratio at 532 nm 

• Color ratio (P 532 /P 1064 ) 

• Depolarization ratio (P perpendicular /P parallel ) 

These lidar measurements are affected by 

• Complex refractive index determined by the aerosol chemical composition 

• Aerosol size distribution 

• Particle shape 

Here we present a study on how lidar measurements of dust aerosol at wavelengths of 532 and 
1064 nm are related to size and complex refractive index. The systematic relationships between 
lidar observables and the dust size and complex refractive index found here that may aid the use of 
space-based or airborne lidars for direct retrieval of dust properties or for the evaluation of 
chemical transport models using forward simulated lidar variables. In addition we show preliminary 
forward simulations of lidar variables based on detailed modeling results 
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Lidar ratio: Scales approximately 
linearly with imaginary part of 
refractive index on log-log scale, 
but slope depends on size. 

Depolarization ratio: Scales 
linearly with real part of 
refractive index for small sizes, 
but dependency is more complex 
for larger particles. 

Color ratio: Scales approximately 
log-linearly with the ratio of 
imaginary part of refractive 
indices at the two observed 
wavelengths, but slope depends 
on size. 


2: Mineral dust optical properties and simulations 

Details of the simulations: 

• Refractive index of minerals and hematite mixtures shown in table 1 

• Refractive index of hematite mixtures calculated with Maxwell-Garnett 
mixing rule 

• Optical properties derived from kernels provided by Dubovik et al. (2006) 

• Assumes mixtures of spheroids with equiprobable aspect ratios between 
0.5 and 2. 

• Lognormal size distributions with 

• r eff values from 0.5-5 pm in 10 steps and 

• v eff values of 0.2, 0.5, 1, 1.5 and 2 

• For minerals with refractive index values outside of the range of the 
optical properties database (m r =1.35-1.59; m^O. 0005-0. 02) extrapolation 
is used 

• Pure hematite is ignored in the simulations because of the extreme 
refractive indices outside the ranges of our database 
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4: Lidar variables dependency on size distribution 


Lidar Ratio vs Effective Radius (Effective Variance = 1.5) 
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Depolarization Ratio vs Effective Radius (Effective Variance = 1.5) 
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Color Ratio vs Effective Radius (Effective Variance = 1.5) 
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Lidar Ratio vs Effective Radius (Effective Variance = 0.5) 
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Depolarization Ratio vs Effective Radius (Effective Variance = 0.5) 
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Color Ratio vs Effective Radius (Effective Variance = 0.5) 


• • llli • • llli-Hema 5.0% 

A A Kaol A A Kaol-Hema 5.0% 

▼ T Smec ▼ T Smec-Hema 5.0% 

♦ ♦ Calc ♦ ♦ Calc-Hema 5.0% 
■ ■ Quar ■ ■ Quar-Hema 5.0% 
► ► Feld ► ► Feld-Hema 5.0% 
4 -4 Gyps 4 4 Gyps-Hema 5.0% 


! f 
t . 


Effective Radius 


Effective Radius 


2 3 

Effective Radius 


Lidar ratio: Increases with size for the highly absorbing minerals, but decrease with size for the weakly absorbing particles. 

Depolarization ratio: For weakly absorbing particles, depolarization ratios increase with size. For absorbing particles this increase with size 
seems to be (partly) compensated by absorption. 

Color ratio: For small particles, color ratios decrease with size, but asymptote to around 1 for weakly absorbing particles, while color ratios 
further decrease with size for absorbing particles. 


Dust model newly developed by Perlwitz et al. 

(ACP 2015, Part 1 and 2) 

• CMIP5 version of NASA GISS Earth 
System ModelE2, 2x2.5 deg., 40 vertical 
layers 

• Prognostic treatment of eight dust types 
in 5 size bins and treatment of iron 
oxides as internal mixtures 

• Includes emission, advection, deposition 
of the minerals, turbulent mixing at 
surface, wet deposition by scavenging in 
and below clouds and condensation 

• Simulations for 2002-2010, using 
prescribed SST and sea ice. Winds 
nudged using NCEP reanalyses . 

• Dust properties are vertically averaged for 
results shown here 

Lidar measurements are simulated for dust- 

only model runs by 

1. Determining appropriate r eff and v eff 
values for each type at each grid box 

2. Determining fraction of area 
contributed by each type at each grid 
box 

3. Use database (panel 3) and area 
fractions to calculate backscatter- 
weighted average of lidar, 
depolarization and color ratios for each 
model grid box (cf. Burton et al., AMT, 
2014). 

Preliminary results (right): 

• Smaller depolarization and color ratios 
and larger lidar ratios near Sahel and 
Australia, reflecting dominance of 
hematite-enriched koalinite and 
smectite, respectively. 

• Relatively low lidar ratios and larger 
depolarization and color ratios of South 
American coast and in China, reflecting 
dominance of feldspar and illite, 
respectively. 
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6: Future work 

• Extend complex refractive index and size ranges in the aerosol optical properties database 

• Simulate vertically resolved lidar variables 

• Calculate lidar variables for model size bins 

• Add modeled non-dust aerosol to global mix 

• Investigate influence of aspect ratio assumptions 

• Include 335 nm wavelengths to investigate possible extra information content 

• Compare with CALIPSO statistics for dusty regions 
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